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ABSTRACT 

We investigate the implications of a bulk rotational component of the halo velocity 
distribution for MACHO mass estimates. We find that for a rotating halo to yield 
a MACHO mass estimate significantly below that of the standard spherical case, its 
microlensing must be highly concentrated close to the Sun. We examine two classes of 
models fitting this criteria: a highly flattened 1/r 2 halo, and a spheroid-like population 
with whose density falls off as 1/r 3 5 . The highly flattened 1/r 2 models can decrease 
the implied average MACHO mass only marginally and the spheroid models not at all. 
Generally, rotational models cannot bring the MACHO mass implied by the current 
microlensing data down to the substellar range. 
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1 INTRODUCTION 



Recent microlensing results from the MACHO collaboration 
suggest that in the context of a spherical isothermal model 
for the Galactic halo, some significant fraction of the halo 
is composed of objects with masses of about 0.4Mq (Al- 
|cock et al. 1996 ). Such masses are consistent with several 
astrophysical candidates for MACHOs - low mass main se- 
quence stars, white dwarfs, and black holes. Each of these 



candidates, however, is subject to a variety of constraints 
that present serious challenges for succesful models in which 
they form a significant fraction of the halo. Low mass main 
sequence stars should be easily visible, and recent direct 
searches for these stars limit their contribution to the halo 



to be less than 3% (Flynn et al. 1996; Graff fc Freese 1996) 



The limits on halo white dwarf stars are much looser. Very 
old white dwarfs would be cool and faint enough to have 
evaded direct detection thus far and yet still be present in 



the halo in significant numbers (Graff et al. 1997 



& Laughlin 1996; dhabrier et al. 1996). The progenitors of 



Adams 



these white dwarfs presumably formed as a very early gener- 
ation of stars. This scenario has its own problems, however. 
The progenitor stars would have produced copious amounts 
of metal enriched gas, which is not seen in the Galactic abun- 



dances (Hegyi et. al. 1986; Gibson & Mould 1997). Ridding 



the Galaxy of this metal enriched gas wo uld require a stron g 
Galactic wind at the appropriate time (Fields et al. 1996). 
In addition, this scenario requires a stellar mass function for 
the halo stars that is peaked at a much higher mass than 



is observed in the disk today. Unless the initial mass func- 
tion is strongly suppressed for low masses we should still 
see the low mass main sequence stars that would have been 
produ ced along with the higher mass white dwa rf progenitor 
stars (Ryu et. al. 1990; Adams & Laughlin 1996). Primordial 
black hole candidates require a fine tuning of the initial dis- 
tribution of density perturbations, although there has been 
some recen t work on black holes formed at the QCD phase 
transition (Jedamzik 1997). 



In this paper we explore an alternative explanation for 
the relatively long term event durations seen by the MA- 
CHO group. MACHO mass estimates are derived from the 
observed event duration, with the assumption of a model for 
the MACHO distribution and velocity structure. The event 
duration is a function not only of the mass of the lens, but 
also of the distance to the lens and its tranverse (with re- 
spect to the line of sight) velocity. Without the additional 
information obtained, for example, by a parallax or binary 
lens event, in which case the velocity of the lensing object 
can sometimes be extracted and an independent estimate of 
the distance to the lens, the mass of an individual event can- 
not be determined directly. The distance along the lensing 
tube and tranverse velocity are known only in a statistical 
sense (for a given model), and therefore only a statistical 
estimate of the mass of a population of lenses can be made. 
However, even given a large number of events, the mass es- 
timate still depends strongly on the assumed distribution 
function (phase-space density) of the lenses. For example, 
MACHOs which are moving through the lensing tube with 
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a transverse velocity that is much less than that in a stan- 
dard isothermal halo model will result in much longer event 
durations for a given mass. 

Brown dwarf stars, with masses below 0.08Mq, repre- 
sent a less troublesome candidate for the observed microlens- 
ing events toward the Large Magellenic Cloud (LMC) but 
appear incompatible with the current data in the contex t 
of a standard isothermal halo model ( |Alcock et al. 1996| ). 
However, the viability of such lens candidates could be very 
different in a model with a non-standard MACHO distribu- 
tion function. 

In order to explore the effects of a non-standard MA- 
CHO distribution function, we allow for a bulk rotational 
component of the MACHO halo and anisotropic velocity 
dispersion. While a large bulk rotation of a spherical halo 
seems unlikely, there is no reason to exclude such a compo- 
nent a priori. Further, the MACHOs, wh ich represent only 
a bout 20% - 40% of the total halo mass (|Gates et al. 1997 



Gates et al. 1996), may be in a more condensed or flattened 



distribution than the CDM halo, and thus it seems reason- 
able to expect some rotational component to their velocity 
profile. 



2 MICROLENSING TOOLS 

In order to compare our models to the MACHO data we 
first calculate the average duration. Note that we are not 
attempting here to determine a detailed model for the halo, 
in which case we would need to consider the event duration 
on an event-by-event basis, but rather to explore the poten- 
tial of non-standard halo distribution functions to affect the 
MACHO mass estimate. The average event duration is di- 
rectly related to the optical depth, rate and MACHO mass 
for a given model. 



imodel — -p , 



In 



Tmodel 



i im y 1M e 

where we have now calculated the rate for solar mass MA- 
CHO's and made the mass dependence explicit. For a given 
model we can then extract the prediction for the average 
MACHO mass, 



(1) 



(2) 
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Calculation of the observed average duration (Appendix 
A) gives a value of 61 days. Due to the low number of e vents 
and the statistics of microlensing ( |Han fc Gould 1995 ), the 
possible error is very large. Nevertheless, for the remainder 
of this paper we will utilize the observed average duration: 
we are interested in how rotation of the halo might effect 
the central value for the MACHO mass estimate. 

The optical depth towards the LMC is evaluated in the 
standard way, 



Rl{l) 



P{1) 



(II 



(4) 



Re{1) = 



AGm l(D 3 - 1) 
c 2 Dl 1 



(5) 



and D s — 50kpc is the distance to the source. 

Calculating the microlensing rate for a moving observer 
and source (the microlensing tube is tumbling) in an ar- 
bitrary MACHO distribution is slightly more complex. For 
each segment dl of the microlensing tube we can calculate 
the corresponding infinitesimal contribution to the rate, dr. 
Let fr{v, I) be the distribution function of the MACHOs in 
the frame of the tube at location I along the tube. For every 
element of the velocity distribution Av x Av y Av z there is a 
contribution to the rate, 



A 3 dF = 2R E (l)dl\v ± \f T (v,l)Av x Av y Av z 
= 2R E {l)dl\v x l\f T {v, l)Av x Av v Av z 



(6) 



where 2RE{l)dl is the cross-section of the segment and I is 
the unit vector pointing along the microlensing tube. Inte- 
grating over all velocities we have 



dV = 2R E (l)dl I \vxl\f T (v,l)d- > v. 

J v 

= 2R E (l)dl < \v ± \ > p(l)/m, 



where < | > is the expectation value of the transverse 
velocity relative to the microlensing tube. Integrating along 
the line of sight we obtain 

r = > I ' i? B (0— < Kl > .//. iM 

m 



It will be useful to consider the average duration due to 
events at a distance I from the Sun, 



to.i(l) 



dr{l) 
dV{l) 
ir Re (I) 



where Re (I) is the Einstein radius at distance I, 



- ., : (9) 

2 < \v± \ > 

where the Einstein radius is evaluated assuming a mass of 
O.IMq- We also consider the fractional event rate, dF/T, as 
a function of the distance I, i.e. the fraction of events that 
occur at any given distance. 



3 EFFECT OF ROTATION 

We start by examining the "standard" halo model, 

2 i 2 

„ „ r + a , . 

P = P0 2 i 2 ( 10 ) 

r z + a z 

where a — 5kpc and ro = 8.0 kpc, with no rotation, flatten- 
ing or anisotropy. Figure 1 shows the transverse velocity with 
respect to the microlensing tube as a function of distance 
from the Sun. The transverse velocity is largest towards the 
ends where the tumbling velocity due to the motion of the 
sun and LMC is greatest, but overall is essentially flat. Fig- 
ure 2 shows the average event duration, again as a function 
of distance, assuming a MACHO halo composed of O.IA/q 
objects. The duration is largest near the center and drops to 
zero for objects close to either end of the microlensing tube. 
The maximum average duration is quite low, only about 43 
days. Regardless of the distribution of event distances, the 
average duration over the entire tube can be no larger than 
this. Scaling to the observed average duration of 61 days we 
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Figure 1. Transverse velocities for the "standard" (solid) and 
rotating (dashed) halos. 



Figure 3. Distribution of the event distances for the "standard" 
halo. 
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Figure 2. Average event durations assuming O.IMq MACHOs 
for the "standard" (solid) and rotating (dashed) halos. 



find a minimum possible average mass of 0.2Mq. The frac- 
tional event rate as a function of distance for this model is 
shown in Figure 3. The actual distribution of events favours 
the first part of the microlensing tube, where the Einstein 
radius is small and durations correspondingly short, leading 
to a much larger estimated MACHO mass, ~ 0.4Mq . Since 
the transverse velocity is relatively flat over the entire dis- 
tance to the LMC, this distribution is driven mainly by the 
distribution of the lenses. 



Next, for comparison, we consider an extreme model 
- an idealized cold co-rotating halo (t> ro t=220 km/s) model 
with zero velocity dispersion (even in the z-direction) . While 
obviously not a realistic model, it does allow us to examine 
the limits on MACHO mass estimate reduction in models 
with a rotational component for the halo velocity structure. 
Figure 1 shows the average transverse velocity of the this 
model as a dashed line. As expected, the transverse velocity 
close to the observer is low. However, note that the trans- 
verse velocity at more distant points is actually quite high: 
the MACHOs sweep past the microlensing tube because the 
LMC velocity is not aligned with the rotation of the Milky 
Way as is the motion of the Sun. The dashed line in Figure 2 
is the average duration of events (for O.IMq MACHOs) for 
the co-rotating halo model. We see a dramatic effect due to 
the reduced velocities close to the Sun: there is a sharp peak 
in the durations of events at very small distances. In gen- 
eral, the durations are increased for distances smaller than 
10 kpc and decreased for distances greater than 10 kpc, re- 
flecting the relative velocities at these distances. There is 
a small region where the durations are above 61 days, the 
observed average duration. This implies that it is possible, 
at least in principle, to obtain average durations matching 
those observed with O.IMq MACHOs, and suggests distri- 
butions which are highly concentrated toward the center of 
the Galaxy as the most likely way to achieve this. 

If a rotating halo is to lead to reduced mass predic- 
tion for the MACHOs, essentially all the lensing must oc- 
cur very close to the observer. The distribution from the 
standard model is too wide and peaks at too large a value 
to fully sample the initial peak in average event duration 
due to rotation. Any benefit from the increased durations at 
short distance will be offset by the decreased durations at 
larger distances. We consider two options for bringing the 
microlensing closer to home: a density fall-off faster than 
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l/r 2 or a flattened halo. Both of these posibilities are quite 
reasonable. The stellar halo population we do know about, 
the spheroid, has a density distribution that falls off roughly 
as l/r 3,5 . The MACHO halo may trace this distribu t ion. On 



the ot her h and, various observations (Oiling 1996; Backett 
et al. 1 )94) suggest that dark halos are flattened, and sim- 



ulations consistently show that the CDM halos generically 
form with departures from spherical symmetry. The dissipa- 
tive baryonic component of the halo, the MACHOs, might 
be expected to be even further from spherical. We know of 
at least one component of the Galaxy, the disk, that became 
significantly flattened before forming stars. In the following 
section we examine these models in detail. 



4 MODELING 

Following the above discussion, we will consider three den- 
sity profiles: 

an isothermal halo model with core radius a 



2 i 2 

r +a 

r A _|_ a A 



a flattened (axis ratio q=0.2) inverse square halo ^| 



P = Po 



2 i 2 

r + a 



R 2 + z 2 /q 2 + a 2 



and a spheroid-like halo 



Po: 



ri 5 + a 3 ' 5 



(11) 



(12) 



(13) 



We fix the core radii at a = 5 kpc for the l/r 2 models 
and a = lkpc for the spheroid model. (We examined the 
effects of varying the core radius for all 3 models and found 
them to be unimportant to our final results.) 

To estimate the MACHO mass implied by these mod- 
els, we need to specify not only the density distribution, but 
also the velocity structure of the halo. We assume that the 
phase space distribution is a simple anisotropic gaussian, 
offset from the origin by the rotation velocity and aligned at 
each point with cylindrical Galactic coordinates.^ We do not 
attempt to construct completely "self-consistent" halo mod- 
els. For the flattened or centrally condensed halos we con- 
sider, such models would need to incorporate the bulge and 
disk, as well as halo self-consistency. This is presently well 
beyond the state of the art given the large uncertainties in 
essentially all Galactic parameters. Furthermore, microlens- 
ing constraints are sufficiently weak (due to low number of 
events) that such sophistication is unwarranted at present. 



* We considered using an Evans model for the halo since these 
come with a consistent prescription for the phase space density. 
However such models cannot produce a halo as flattened as we 
wished to explore, at most obtaining an equivalent flattening in 
the density of R3 1/3. 

t Most of the lensing in our models takes place close to the Solar 
position, where a velocity ellipsoid aligned with Galactic cylin- 
drical coordinates is a reasonable approximation, and over a rel- 
atively narrow range in distance where the anisotropic velocity 
dispersions should not change radically. Thus, while these are not 
the most general set of assumptions, they are sufficient for our 
purposes and greatly simplify the calculations. 



The Jeans equations for the above time-independent 
distributions simplify to, 



9(pv%) , P 



9$ 



dR 



djpvl) d$ _ Q 
dz dz ' 
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(15) 



where <!> is the potential of the total halo. We further assume 

that ~ and ^j- « at points in the halo of interest 
for this work. Although this is probably a poor assumption 
for detailed modeling, we accept it in order to focus on the 
effects of rotation on the microlensing mass estimates. We 
then have 



d In p \ , 2 . — 2 

- 1 ) <TR + °~4> + v<t> 



d\nR 



= R 



where a 



R 



(£)•* - - 

a 2 = v 2 and a% + 



9£ 
dR 

fir 

OZ 



(16) 

(breaking it up 



into random and bulk flows). The derivative ^ is just the 
R-component of the force at R. Since we are roughly at the 
equator, i?ff| is simply the square of the circular velocity 
at R, v 2 (R) ps (220km/s) 2 

To complete this set of equations we need to find §§. We 
assume that the bulk of the halo, and thus the major contri- 
bution to the potential, comes from the non-baryonic com- 
ponent, whose distribution is assumed to be that of equation 
hi with core radius onb = and arbitrary flattening qnb in- 
dependent of the MACHO halo. It can be shown (Appendix 
B) that to first order in z 



9$ 

dz 



Inb 



z 

R 2 



(17) 



near the equatorial plane. This gives 



2 2 

o- z = ~pv c 



Inb 



sin 1 y/l - qg, B 



The lowest value for a 2 is obtained for a spherical (q=1.0) 
non-baryonic halo. A larger value of the vertical velocity 
dispersion will drive the transverse velocities up, the event 
durations down, and ultimately yield a higher mass esti- 
mate. Since we are interested here in models which lower 
the MACHO mass estimates we take qnb = 1 and find 



\ d\nR J 



1 I o\ + a 2 



(19 



(19) 



-PVo^ 



As we will see, the precise form of the relationship be- 
tween the rotation and velocity dispersion is not important 
for our conclusions. We use equation as a general guide 
to which combinations of rotation and velocity dispersions 
are reasonable. 



5 RESULTS AND DISCUSSION 

We plot our results for the various halo models in figures 
4 through 6. The solid contours correspond to the pre- 
dicted MACHO mass estimate, corresponding to the aver- 
age MACHO event duration of 61 days, as a function of the 
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Figure 4. Average MACHO mass for an isothermal halo based Figure 5. Same as Figure 4 for the flattened 1/r 2 halo. Countours 

on the MACHO collaboration average event duration. Countours are 0.15, 0.2, 0.25, 0.3, ..Mq from left side, 

are 0.15, 0.2, 0.25, 0.3, ..Mq from left side. 



(one-dimensional) velocity dispersion [] and rotation speed. 
For reference, the standard non-rotating halo with isotropic 
maxwellian velocity dispersion has a one-dimensional dis- 
persion a — 156km/s. 

The models within ±10% of our a — v Iot curve (equation 
|is| ) lie in the shaded region between the dotted lines. Models 
below this region are unlikely to have sufficient support. The 
trade-off of dispersion velocity for rotation velocity can be 
seen clearly. 

Figures 4, 5 and 6 show the predicted MACHO masses 
for the spherical 1/r 2 , flattened 1/r 2 and spherical 1/r 3 ' 5 
halo models. In no case does the predicted mass for the 
current MACHO event duration go below O.25M0 in the 
allowed region indicated by equation [H^ . As expected fol- 
lowing the discussion in section 3, the spherical 1/r 2 model 
actually predicts increasing MACHO mass estimates as it 
becomes more rotation supported. Rotation increases the 
velocities along the microlensing tube far from the Sun lead- 
ing to shorter event durations for a given MACHO mass, or 
conversely, a larger MACHO mass estimate for a given ob- 
served event duration. Most microlensing in the standard 
halo model occurs at distances where the velocities are in- 
creased when a rotational velocity component is added and 
hence the mass estimates for a given event duration increase. 

In contrast, in highly flattened 1/r 2 halo models rota- 
tion does lead to reduced mass estimates. However, the im- 
provement is not as dramatic as might be hoped. The reason 
is a little subtle. In the absence of rotation, the microlensing 
is concentrated in the first part of the tube by the flattening 

t In the discussion that follows we take a = an = a^. We have 
also considered models with a constant anisotropy an = ota^,, 
where 0.3 < a < 3.0; our results are essentially unchanged for 
these models. 



of the halo. We can see that predicted mass estimates for 
the flattened halo model with no rotation are a little higher 
than for the spherical model since in this region the du- 
rations for the non-rotating model are short. However, the 
fractional rate is also dependent on the average transverse 
velocity of MACHOs. Because a flattened halo has a lower 
z-velocity dispersion, when rotation is added the transverse 
velocity can become very low close to the observer. This is 
precisely why rotation might be expected to lower the pre- 
dicted MACHO mass estimates. However, the low transverse 
velocity also suppresses the event rate close to the observer 
and shifts the distribution of events to larger distances where 
the transverse velocities are higher. These two effects com- 
pete, with the end result that the reduction to the MACHO 
mass estimate due to rotation is only modest. 

Our second hope for taking advantage of rotation, a 
more concentrated 1/r 3 ' 5 spheroid, is also disappointing. 
Supporting a spheroidal halo rotationally leads to very lit- 
tle change in the MACHO mass estimate. Rotation fails to 
reduce the mass estimates in the spheroidal case for reasons 
slightly different than those for the flattened halo. For a 
spheroidal MACHO distribution, transverse velocities never 
get very low, even close to the observer. Hence, little benefit 
is derived from the central concentration. Such a spheroidal 
model also has serious difficulties producing a high optical 
depth without ruining the flatness of the Galactic rotation 
curve. 

In this section we have explored three general classes 
of halo models, which have been unable, even with the ad- 
ditional of a generous rotational component, to reduce the 
MACHO mass estimates for the current observed event du- 
rations to masses in the substellar regime. The failure of ro- 
tating halos to significantly reduce MACHO mass estimates 
can be summed up as follows. First, the virial theorem im- 
plies that a typical velocity for a particle in an extended halo 
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Figure 6. Same as Figure 4 for the spheroidal halo. Countours 
are 0.2, p.25, 0.3, ..M® from left side. 



will be roughly the same regardless of whether the halo is 
or is not rotationally supported. Velocities cannot be made 
arbitrarily low. Secondly, the direction to the LMC is "awk- 
ward", that is, it respects none of the symmetries of the 
Galaxy. This has several implications for attempts to lower 
the MACHO mass estimates. Bulk velocities due to rotation 
cannot be used to cancel out some of the motion of the MA- 
CHOS through the lensing tube except possibly very close to 
the observer. An effort to concentrate the MACHOs in this 
region, however, is self-defeating in the models we explored. 
The event rate in such models is dominated by MACHOs 
moving through the tube further away from the observer, 
and the net effect on the MACHO mass estimate is not suf- 
ficiently large to reach the substellar regime with current 
data. Further, the direction of the LMC makes it impossi- 
ble to arrange, in any natural way, an anisotropic velocity 
dispersion such that the net t ransverse velocity through the 
tub e is significantly reduced (Gyuk & Gates in preparation 
19970 



6 CONCLUSIONS 

Scenarios in which brown dwarfs populate a standard spher- 
ical, non-rotating halo in significant numbers predict event 
durations significantly (more than a factor of 2) shorter than 
observed in the current MACHO data. We have explored 
halo models with a bulk rotational component in an attempt 
to lower the predicted mass estimate for the MACHO events 
to the substellar regime. We find that unless essentially all 
of the lensing takes place within 1-2 kpc of the Sun and 
residual velocity dispersions are very low, < 30km/s, it is 
not possible to reduce the MACHO mass estimate in these 
models to the brown dwarf mass range. Such a configuration 
no longer resembles a halo distribution for the MACHOs. 
A highly flattened rotating halo was our most successful 



model in reducing MACHO mass estimates, although even 
this model did not reduce the microlensing mass estimate 
below about 0.25M®. This suggests that a thick disk con- 
figuration with a steeper (exponential) density fall off away 
from the Galactic plane might be able to reduce the mass 
estimates further. A high rotation velocity is also more nat- 
ural for such a distribution than for a more extended halo. 
Preliminary work shows that a thick disk can indeed reduce 
the mass estimate, but only to around 0.15 — Q.2Mq. We 
will report on this in more detail in a forthcoming paper. 

In summary, the prospect for a brown dwarf microlens- 
ing halo appears dim given the current data, and the ques- 
tion of what and where the MACHO lenses are remains 
unanswered. 
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APPENDIX A : CALCULATION OF AVERAGE 
EVENT DURATION 

From the observed events we can construct the observed dis- 
tribution function of time scales, v(t). The true distribution 
will then be p(i) = u(t)/e(t)/ J v(i)/e(i)di. We want to cal- 
culate 



(20) 
(21) 



t = J tp(t)dt 

J ii/(t)/e(t)dt 
J v{i)/e(i)di 



which requires knowing the function v(t). Our best guess for 
this will be to let 



i 

Substituting in gives 



(22) 
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(23) 



Using the events and efficiencies towards the LMC ( Alcock 
et al. l!)96) we obtain t — 61 days. 



APPENDIX B : VERTICAL POTENTIAL 



Following (Binney & Tremaine 1987), the z-component of 
the gravitational force generated by a density distribution 
p — p(m 2 ) where m 2 — R 2 + z 2 /q 2 , is given by 

-g^wr^r j im2)zdm — , (24) 

9 * Jo m[S-(l-q 2 )} 3/2 

where we keep only the leading order term in z, in order to 
obtain a lower limit on the estimate for a z . 

Assuming a density p — pom^/rn 2 we then arrive at 

oz R- 1 
Rewriting this in terms of the circular velocity and flattening 
we find 

" < " (26) 



dz y/q 2 /{l-q 2 )suT 1 {l-q 2 )K 1 ' 
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